The impact of the spatial legacies of spruce budworm defoliation on wildfire burn severity and post-fire recovery

Background

Boreal forests are highly complex ecosystems that are shaped by the interaction between forest structure and composition, topography, climate and weather, and disturbance regimes. This complexity causes general uncertainty around what are the particular drivers of ecosystem processes, such as wildfires. Wildfires are the primary disturbance agent that allows for forest regeneration and helps to maintain forest structure and function. Because of the central role that wildfires play in the boreal forest, it is therefore crucial to understand where they occur, why they occur and what factors influence their behaviour. Many researchers have posited that insect disturbance such as spruce budworm defoliation, which is another widespread disturbance agent across the boreal forest, has a strong impact on wildfire likelihood and behaviour by changing the abundance, structure and connectivity of combustible fuels. However, wildfires are relatively uncommon in the boreal forest, occurring on average between 50 and 200 years and the spatial legacies of insect disturbance persist just for a short window (~15 years). Additionally, both wildfires and spruce budworm defoliation affect a vast amount of area annually (give stats), therefore it is difficult to detect a small signal increas noise….

To date, current evidence on the impact of the spatial legacies of defoliation is mixed, reporting positive, negative and neutral effects. This uncertainty is in part driven by the ability to detect where particular changes in fuel or legacies of defoliation have overlapped with wildfires. Recent advances in remote sensing offers exciting prospects to utilize novel spatial data combined with thoughtful experimental design to reduce uncertainty around the effects of spruce budworm defoliation on fuels and wildfire behaviour. We are now able to remove or reduce uncertainty through better experimental design and approaches to big data problems. Here, we used a paired design within fires that have both areas that were affected by insect defoliation and not affected to attempt to resolve this uncertainty,

Questions

The overarching question is: how do the spatial legacies of spruce budworm defoliation influence wildfire burn severity and forest recovery in the boreal forest?

Specifically we will address the following questions:

  1. How does spruce budworm defoliation influence burn severity?
    1. Is there a difference in median severity between defoliated and non-defoliated fires?
    2. Is there a difference in severity extremes between defoliated and non-defoliated fires
    3. Is there a difference in variability in burn severity between defoliated and non-defoliated fires
  2. How does spruce budworm defoliation influence post-fire recovery?
    1. Is there a difference in forest recovery 10 years post-fire between defoliated and non-defoliated fires?
    2. Is there a difference in forest recovery between 1-5 years and 6-10 years between defoliated and non-defoliated fires?
  3. Is the strength of the relationship between defoliation intensity (synthetic index) and burn severity or forest recovery greater than independent metrics of defoliation?
  4. Does defoliation intensity (synthetic index) capture more variation in burn severity and forest recovery relative to defoliation metrics modeled independently?

Methods

This study focuses on wildfire that burned in the boreal shield ecoregion of Ontario. We restricted our analysis to wildfires that burned between 1986 and 2013. 1985 was chosen as the first fire year in our dataset because it is the first year that remotely sensed burn severity data are available. Similarly, our analysis of forest recovery required landsat imagery for every year for 10 years post-fire, therefore 2013 is the last year where we were able to calculate remotely sensed imagery. We identified wildfires that were defoliated up to 15 years prior to the year of fire. We chose to limit the defoliation to within 15 years because previous studies have shown that beyond 15 years, defoliation is not significant due to enhanced regrowth. We only selected fires that contained both defoliated and non-defoliated areas. Once identified, we split each fire into two distinct areas. This resulted in a total of 34 (n = 68) fires containing both defoliated and non-defoliated areas.

Structural Equation Model

The second part of the analysis looks at the effects of defoliation intensity on burn severity and forest recovery following fire. Defoliation alters fuel structure and fuel load and is expected to cause influence fire behaviour. Previous studies in Ontario have shown that defoliation can alter fire behaviour for up to 15 following defoliation [@james2017], this period is commonly referred to as the “window of opportunity”.

Spruce budworm defoliation is multivariate in nature. The effects of spruce budworm defoliation can be described by the presence and absence of defoliation, the time since defoliation to the time of fire and the number of years that an area was defoliated (a proxy for defoliation load). It is common practice for studies that concerned with the effects of spruce budworm defoliation on wildfire to use a since predictor, or build multiple sequential models to incorporate the multivariate nature of defoliation [@james2017 @james2011 @candau2018].

Structural equation models allow for the multivariate nature of defoliation by creating composite variable that allows us to combine the variation explained by multiple predictors into a synthetic index. Using structural equation models, we will combine create a synthetic index called defoliation intensity which combines defoliation presence/absence, time since defoliation and number of years defoliated.

We will fit these models to data of consisting of 34 paired defoliated/non-defoliated fires that burned in Ontario between 1986-2012. We will then model the causal effects of defoliation intensity on burn severity and forest recovery following fire. Forest recovery is measured by taking the mean composite of the normalized burn ratio for a given summer for 10 years for following fire. We then use a Theil-Sen regression to calculate the sen’s slope across the 10 years which we call the slope of recovery or forest recovery.

In our model, defoliation intensity will effect for burn severity and forest recovery. Below is the directed acyclic graph (DAG) of our a priori causal hypothesis.

In this DAG we see that time since defoliation, and number of years defoliated are combined to form the synthetic index (defoliation intensity) which is then used to predict the effects on burn severity and forest recovery. Defoliation intensity will indirectly effect forest recovery by altering stand structure and forest composition prior to fire, which will impact the number and distribution species and their seed banks that are able to establish post fire, thus affecting the recovery trajectory. Similarly, defoliation intensity will effect burn severity by either increasing or decreasing severity which will lead to a change in the recovery trajectory by either increasing or decreasing the slope of recovery.

Preliminary results

How does spruce budworm defoliation influence burn severity?

## Loading required package: gt
Results from Pairwise T-test on severity between defoliated and non-deoliated
Test t-value df P-value
Median Severity -12.1647870670551 65 <0.0001
Extreme Severity -25.3390680767449 65 <0.0001
Variability in Severity -2.83043157076872 65 0.006

Severity visualization #### How does spruce budworm defoliation influence post-fire recovery?

Results from Pairwise T-test on Slope of Recovery Between Defoliated and Non-Deoliated
Test t-value df P-value
10yr post-fire -7.22459615778699 65 <0.0001
1-5yr post-fire -3.09145675465788 65 0.003
6-10yr post-fire -3.43414952796882 65 0.001

Slope visualization

targets::tar_read(lolliplot_slope)

Is there a difference in the direction of the effects ?

chisq_table(targets::tar_read(chisq_result))
Significant Difference in Direction of Trend between Paired Fires
Results from chi-square tst
Independent Variable of Interest df P-value
Median Severity 1.515 1 0.218
Extreme Severity 0.545 1 0.460
Varibility in Severity 7.333 1 0.007
10 yr slope of recovery 0.061 1 0.806
1-5 yr slope of recovery 2.970 1 0.085
6-10 yr slope of recovery 0.545 1 0.460

Does a synthetic index of defoliation intensity capture more variation in burn severity and forest recovery relative to defoliation metrics modelled using an interaction? Are the strengths of the effects different?

First, we fit repeated measures linear mixed effects models with an interaction terms between time since defoliation and cumulative years of defoliation for each response variable:

\(Y_{i} = \beta_0 + \beta_1X_{1i} + \beta_2X_{2i} +\beta_3X_{1i}X_{2i} + u_{i} \epsilon_{i}\).

Where Y = the response variable of interest and i = for the i th fire, \(\beta_1X_{1i}\) = time since disturbance (tsd), \(\beta_2X_{2i}\) = cumulative years of defoliation (cyd), \(+\beta_3X_{1i}X_{2i}\) = interaction between tsd and cyd, \(u_{i}\) = random intercept for paired fire id.

Median Severity

## Warning: package 'dplyr' was built under R version 4.3.1
Median Burn Severity
Results from repeated measures hierachical linear mixed model fit in nlme
term estimate std.error df statistic p.value
(Intercept) 137.918992 14.845430 32 9.2903331 0.000
tsd -6.008260 3.521365 30 -1.7062304 0.098
cumltve_yrs -1.271297 7.934572 30 -0.1602225 0.874
tsd:cumltve_yrs 2.757749 1.519297 30 1.8151480 0.080

Extreme Severity

Burn Severity Extremes
Results from repeated measures hierachical linear mixed model fit in nlme
term estimate std.error df statistic p.value
(Intercept) 416.938629 20.057456 32 20.7872144 0.000
tsd -12.824795 3.351523 30 -3.8265572 0.001
cumltve_yrs -3.092247 7.563267 30 -0.4088507 0.686
tsd:cumltve_yrs 3.497500 1.495899 30 2.3380595 0.026

Variability in Severity

Variability in Burn Severity
Results from repeated measures hierachical linear mixed model fit in nlme
term estimate std.error df statistic p.value
(Intercept) 0.92211373 0.098091067 32 9.4005881 0.000
tsd -0.07496247 0.021369826 30 -3.5078650 0.001
cumltve_yrs 0.01831034 0.048182540 30 0.3800203 0.707
tsd:cumltve_yrs 0.01222242 0.009355431 30 1.3064518 0.201

Slope of recovery 10 year

Slope of Recovery 1-10 Years Post-Fire
Results from repeated measures hierachical linear mixed model fit in nlme
term estimate std.error df statistic p.value
(Intercept) 4.75171721 0.82529339 32 5.7576097 0.000
tsd -0.03680817 0.08095987 30 -0.4546471 0.653
cumltve_yrs -0.05826139 0.18279534 30 -0.3187247 0.752
tsd:cumltve_yrs 0.01962360 0.03654637 30 0.5369508 0.595

Slope of recovery 1-5 year

Slope of Recovery 1-5 Years Post-Fire
Results from repeated measures hierachical linear mixed model fit in nlme
term estimate std.error df statistic p.value
(Intercept) 4.69897749 1.8953407 32 2.4792258 0.019
tsd -0.27094399 0.2104410 30 -1.2875057 0.208
cumltve_yrs 0.23714749 0.4751126 30 0.4991396 0.621
tsd:cumltve_yrs 0.05889672 0.0948599 30 0.6208810 0.539

Slope of recovery 6-10 year

Slope of Recovery 6-10 Years Post-Fire
Results from repeated measures hierachical linear mixed model fit in nlme
term estimate std.error df statistic p.value
(Intercept) 5.77401694 2.27665532 32 2.5361841 0.016
tsd 0.09521720 0.17509967 30 0.5437886 0.591
cumltve_yrs -0.29121359 0.39538353 30 -0.7365344 0.467
tsd:cumltve_yrs 0.03183047 0.07918944 30 0.4019534 0.691

Significance

Land managers in the boreal forest region require information about burn severity to address immediate hazards and outcomes resulting from wildfires. Understanding how previous disturbance history shapes patterns of burn severity is needed to guide management activities and ensure downstream ecosystem provisioning when threatened by multiple disturbances. Results from this chapter will provide a framework for modelling the effects of spruce budworm defoliation on wildfire severity and help managers identify areas that are at the highest risk of losing valuable ecosystem services.